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Unoer  this  grant  the  researchers  have  Drimarily  Investigated  two  classes  of  e letromagnet  lc 
c-otileins.  The  first  contains  the  ouant i tat Ive  desci'lDtion  of  nicrowava  heating  of  dlsoersive  and 
conductive  materials.  Such  oroblema  arise,  for  example,  when  biological  tissue  are  exoosed. 
accidentally  or  purposefully,  to  microwave  rad'ation.  Other  instances  occur  in  ceramic 
Dr-jcess’ng.  such  as  sintering  ang  microwave  assisted  chemical  vapor  infiltration  and  other 
industrial  drying  processes,  such  as  the  curing  of  paints  and  concrete.  The  second  class 
characterises  :hc  scattering  of  microwaves  by  complex  targets  which  possess  two  or  more  disparate 
length  and/or  time  scales.  Spatially  complex  scatters  arise  in  a  variety  of  ape  1 i cat  ions ,  such 
as.  large  gratings  and  slowly  chariging  gviidins  structures.  The  former  arn  useful  in  developing 
microstrip  energy  couplers  while  the  later  can  be  used  to  model  anatoin*cal  subsystems  (e  g.  the 
open  uuidinvi  structure  composed  of  two  lags  and  the  adjoining  lower  torsol  Temuorally  comolex 
targets  occur  in  applications  involviiig  dispersive  media  whose  relaxation  times  differ  by  orders 
of  magnitude  from  thermal  and/or  electromagnetic  time  scales.  I'or  betn  cases  tne  mathematics' 
description  of  the  problems,  give  r’se  to  comp  1  iceted  i  1 1 -cend  1 1  lored  boundary  value  orcc'ofis.  i 
whose  accurate  solutions  recuire  a  blend  of  both  asymptotic  technlaues.  such  -js  nultiscale  ntthode  1 
and  matched  asymptotic  expansicns.  and  numerical  methods  t  ncorpora  i  ■.  ng  racliat'.cn  bouraarv  , 
conditions,  such  as  finite  differences  and  finite  oloments.  _ 
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Final  Rviport 

Contraci  Number  AFOSR  91-0252 
y  pplied  Mathematical  Problems  in  Modern  Electromagnetics 

Prepared  for: 

Air  Force  OfUce  of  Scientific  Research 


Research  Activities 


UncUr  tlu'  ciirrcnt  grant.  \v<'  hav«‘  primarily  inv<‘stigat<'<l  two  cIjusss'k  of  <'loc- 
troinagui't ic  j)rol)li'ms.  Tlu'  first  <'ontaius  the  <jtiantitativ<-  (l(‘scrii>tiou  of  mi- 
l•rowa\•^•  h(‘:itingof  disporsira' and  rondnctivc  materials.  Sticli  prnhlnus  arise, 
lor  I'xample.  wIkmi  hiologietil  tissue  iire  expo.s«-d.  in‘<-identally  or  ptirposefully, 
to  minowiive  raditition.  Other  iiistime<‘s  ocrnr  in  ceramie  i)ro(essing.  smh 
as  sinteiim;  an<l  mi<'idw:ive  assisteil  chemifid  vapor  iniiltiiition  and  otlau 
iialiistiial  drying  juoeesses.  sni’li  it.-'  the  I'uriiig  of  piiints  and  «a)nerete. 

The  '•(■foiul  elas>  characttui. the  seattering  »)f  mirrowtives  hy  ;-omj)l<’X 
taiget-.  whii'h  p<issev>  two  or  moti-  dispaiate  length  ancl/oi  time  scales.  Sp:i 
tially  i-oni]>h  X  ^eatteieis  arise  in  a  vaiiefy  of  applications,  sueli  iis.  large 
"ratings  and  slowly  .■•iianging  guiding  structrires.  The  former  arc  useful  in 
develiiprng  rnrerosttrp  I'liergy  ••ou|ders  while  the  later  can  he  used  to  rirmlel 
anaiomti  al  'uh-\ sterrrs  le  g  the  o|M-tt  guiding  structure  cornpo.sed  of  two 
!eg>.  and  the  adjoruitig  lower  torso),  lem]>otally  complex  targets  oei-irr  in  ap 
pin  atrori'  rn\ohitig  dispersivi-  media  whose  relaxation  times  dillei  hy  orders 
ol  uiagnttinle  liom  thermal  and/or  elect iomagn<‘tic  time  scales. 

[■or  hot  h  lases  the  nial  heinat  ical  <lesciiption  of  thi'  piohlems  gi\-es  i  ise 
to  corn] ilicai I  < I  ill  I'oiidit lotted  houmhity  Vidtie  piohlems.  whose  accurate  .so 
lutioiis  i((|uiii  a  tileiid  of  hotli  :i.symptotic  techni(|ue.s.  such  as  multiserde 
lin  t  hod  ■  and  mateln  d  .'isyin|iiotic  expansions,  and  iiumei  ical  met  hods  incoi- 
p<ti tiling  radiation  houinhuy  conditions,  su<-h  as  linite  diii’ert'iiees  and  Unite 
elnilent  s 

W  <  tilso  hegaii  some  pielimiiiiiiy  Work  on  comptiting  tin'  shapi'  of  a  thrc'e- 
diinensional  target  from  i iitie  dependent  scatti'iing  datii.  I'he  aiiplicatirius  of 
iineiM'  inohlems  of  this  type  are  numerous  in  both  indu.st.riiil  anti  biological 

applications. 

l-'inallv.  wi  hegan  to  develop  eificient  numerical  trlgoritlmis  for  determining 
how  ineideni  jiulses  scatter  from  a  variety  of  targets. 

Hefoie  sunini.ii'i/.ing  our  ju-ogress  we  note  tlurt  four  Apix’udice.s  ar(^  in¬ 
cluded  with  this  re])ort.  Appi'udix  A  contains  n'ferenct's  that  are  portiiient 
to  out  research  summari/.ed  helow  and  are  noted  in  the  text  by  Roman 
.N’liinetals.  .\ppendix  B  eontains  a  bilrliography  of  paper.s  published  or  sub¬ 
mitted  lot  ptihlieatiou  that  have  been  sui)ported  l;y  this  contract.  Appendix 
('  eontains  a  list  of  pii'sentations  that  have  been  made  during  the  contract 
pi'iiod.  .And  .Ajipendix  D  contains  a  list  of  recognitions  obtained  and  other 
a  ■cuiiiplisliineiits  made  dming  this  period. 


1,  optimal  Power  Absorption  in  Lossy  Targets. 

VVluni  n  tiuK'  haiiuouic  elcctroiiiagiictic  wave  impinges  upon  a  lossy  target 
a  ])()itioii  of  the  incident  wave  is  absorbed.  This  absoiption  is  due  to  ohmic 
hc'ating  when  the  mat('rial  is  nu'talli<-,  un<l  i.s  duo  to  dipolar  heating  when 
the  material  is  a  dieleotiie.  In  either  ease  tin'  power  absorlx'd  by  the  target 
is  givj'ii  by 

(1  )  ^  ^  Ji 

where  fT  is  the  effective  coiahict ivitv,  E  is  the  «'lectrie  field  within  the  tar¬ 
get,  and  D  is  the  compact  region  oecupieil  Ity  the  scatterer.  An  important 
•  luestiou  to  l)iotuedical  res<'iirehei s  interested  in  h<’ating  lisstu'  and  indnstriiil 
lesearehers  inteic-sted  in  material  piocessing  is  thi'  exisleinc  iin<l  comitnta 
tioii  of  iin  optimal  conductivity  where  tlie  material  heats  most  efficiently. 

Koi  sniiill  i7  it  is  appjireiit  from  ( 1)  an<l  Ma.xwell's  <‘«juations  tliat 

(2)  f  lE„|^/.rWvd: 

-  Ji) 

where  Ed  is  the  electric  tieM  in  tlw  tauget  when  n  —  ().  Tb.us,  Q  approiiehes 
y,<'ro  lineiuly  with  cr.  On  the  other  hiind  its  fr  —*  s.,  the  elc'clric  field  becomes 
exponent iailly  small  outside  of  a  very  thin  botmdary  layer  (a  few  skin  depths 
in  thickiK'ss)  so  tlnit  tlie  iib.sotbed  po\v(‘r  given  by  (1)  iip))roaches  zc'ro  in 
this  limit  too  [i-iij.  Since  the  fiiuetions  :ipp«'iiring  in  (1)  an*  different iiil  with 
respc'Ct  to  a.  ii  imiximuin  absorbed  power  Q*  «>ecins  for  iin  optimal  rr  —  a* . 

The  functioiiiil  dependeia'c  of  <7*  iind  Q*  upon  tlu'  various  electrical  iind 
geometricid  piU'iimeti'is  of  the  tanget  will  lx-  coiuplieaited  and,  in  generaih  aic- 
ce.ssiblc  only  through  I'xtensive  numerieail  scaittering  coinptitaitio’us.  In  Ref- 
{’leiice  1  we  Inwe  pi'e.s<-nted  a  si\ni>l<'  aipjaroximate  method  whii'h  explicitly 
gives  this  relaitionslup  auid  signifieaintly  reduces  the  nuinlaei  of  computavtions. 
We  havve  ehosen  the  problem  of  TM  .scaittering  from  ai  finite,  hcinogeneous, 
and  1(  sy  dielectric  shil)  to  demonstrate  tlu'  method,  because'  we  can  compare 
the  exact  results  with  our  apiaroximaitions. 

Specifically,  we  Inive  obtaiiiu'd  huge  cr  an<l  smavll  a  (?x])ainsions  of  the  solu¬ 
tions  directly  from  Maixwell's  ecpuitious  aind  naive  computed  the  correspond¬ 
ing  limits  of  (1).  We  have  used  a  pade  [iii]  aiiiproximaition  to  blend  these 
results.  This  aipproximation  Inus  both  limits  built  into  it  as  well  ais  ai  smooth 
transition  between  them.  The  validity,  aiceuraiey,  and  refinements  of  this  aip- 
[jroaich,  as  well  ais  its  extension  to  higlier  dimensions,  au'('  under  investigation. 
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We  have  recently  begun  the  study  of  an  analogous  problem  in  which  a 
Debye  medium  absorbs  microwave  energy.  The  idea  is  to  blend  the  low  fre¬ 
quency  and  high  frequency  expansions  of  the  electric  field  together,  using 
padc  [iii]  and  Galerkin  [iv]  approximations,  to  obtain  an  approximation  that 
is  uniform  over  a  large  bandwidth  of  frequencies.  Preliminary  calculations 
on  the  j^ade  approximation  are  quite  encouraging  and  refinements  of  these 
approaches,  as  well  as  its  extension  to  higher  dimensions,  are  under  investi¬ 
gation. 

2.  Microwave  Heating  of  Ceramic  Materials. 

We  have  modeled  and  analyzed  the  heating  of  a  ceramic  slab  [2,4]  under 
TM-ilhuniiiation  in  the  small  Biot  nmnbcr  regime.  The  temperature  dis- 
tril)ution  is  almost  si)atially  uniform  in  this  limit  and  its  evolution  in  time 
is  gov('rned  by  a  first  order  nonliiuxir  amplitrule  equation.  This  equation 
admits  a  time  indei)endent  solution  which  is  a  multivalued  hmction  of  the 
inier()wav(‘  powc'r.  Tlu'  graph  of  this  steadj'  temperature  as  a  function  of  the 
microwave  power  gives  an  S-shaped  response  curve  when  the  electrical  con- 
duetivity  is  mod<’led  ('ither  as  an  exponential  function  of  t('inpcrature  or  an 
Arrlu'iiius  law.  The  dynamics  of  the  heating  process  arc*  deduced,  from  the 
amplitude  equation  and  the  multi vahu'd  resi)ons(‘,  and  itre  dependent  upon 
the  microwave  power  and  initial  conditions.  For  certain  initial  conditions  and 
power  levels  the  system  evolv(\s  to  the  upper  hruneh  of  the  resjunise  curve 
which  corn'sponds  to  tln'iiual  runaway.  Other  initial  conditions  and  power 
h'Vcls  force  th('  system  to  evolve  to  the  lower  hraneh  and  a  sah'  sint(-ring 
t<'mp<'rature. 

This  heating  process  can  be  controlled  in  .some  senst'  by  varying  the  uii- 
erowave  i)ow<-r  in  time  at  a  rate  comnK'usuratc  with  the  tluninal  changes 
(v).  Specifically,  tlie  i)o\ver  is  allow<'d  to  chang«‘d  in  an  expoiu'ntial  fashion 
from  a  higla'i  to  a  lower  power  level.  This  n'lntionsluj)  is  fiiriu'd  into  a  dif- 
fc'iential  (-([uatiou  wliicli  is  appended  to  the  aini)litude  (’(piation  to  form  an 
autonomous  system  of  tlu'  second  onh’r.  Thi.s  system  is  analyzi'd  using  a 
pluise- plane  methods.  The  analysis  shows  the  existema-  of  a  sttibh'  manifold 
wdiich  divides  th<'  pliii,s<'  phnie  into  two  pju  t.s;  Tiitjectoi  i('s  in  the  region  itbovt' 
thi.s  curve  correspond  to  rtnniway  heating  while  thos('  beh'w  yi<‘ld  sttibh'  sin¬ 
tering.  Various  heating  scemuios  are  presented  and  iiiscu,ss<'d  in  Iti’fi'reiict' 
2. 

We  hav('  considered  a  slightly  more  eonipli<'at('d  control  procc'ss  [G]  in  which 
the  r;ite  of  chtinge  of  the  micianvave  power  in  time  is  etpuil  to  the  sum  of 
two  terms;  the  first  is  itroitortional  to  the  po\v«'r  :uid  the  second  to  the 
t('mpcr:tturc.  When  th('  second  term  i.s  ncgh'cted  tlu'  i)ower  is  an  expoinnitial 
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function  of  time  that  increases  from  a  given  initial  value  to  a  final  value  which 
is  chosen  to  be  consistent  with  the  sintering  temperature  [2,4]-  For  a  wide 
value  of  parameters  the  analysis  shows  that  the  system  docs  evolve  to  the 
sintering  temperature.  Unfortunately,  the  trajectory  of  the  system  passes 
through  temperatures  which  are  above  the  melting  point  of  the  material  and 
this  fact  suggests,  in  the  context  of  the  one-dimensional  model,  that  the 
sample  is  destroyed.  However,  when  both  terms  are  present  in  the  control 
equation  and  they  are  properly  chosen,  the  analysis  shows  that  the  sintering 
temperature  can  be  reached  without  sucli  an  overshoot.  Improper  choices 
give  rise  to  temperature  overshoots  and  even  relaxation  oscillations. 

3.  Microwave  Heating  of  Dispersive  Media. 

We  have  recently  extended  arc  one-dimensional  work  [8]  to  study  the  heat¬ 
ing  of  a  compact  dispersive  target  by  a  pulsed,  plane  microwave.  The  disper¬ 
sive  character  of  the  medium  is  described  by  a  Debye  model  and  the  conduc¬ 
tive  naturt!  is  modeled  by  an  ionic  conductivity.  The  electrical  and  thermal 
paranietcrs  are  allowed  to  depend  upon  temperature  which  gives  rise  to  a 
highly  nonlinc'ar  initial,  boundary  value  problem.  The  governing  equations 
are  av<'ragod  in  the  limit  as  uT/i  oo  wlu're  T//  is  a  characteristic  tinu'  at 
which  heat  diffus<'s  in  tlu*  targc't  an<l  u;  is  th('  microwavi'  carric'r  ir('quency.  A 
two-step  algorithm  is  proposed  for  the  munerical  solution  of  tlu'se  eq\iatious. 
When  convection  is  w(!ak,  th<'  algorithm  converges  very  slowly.  Howevtu', 
this  problem  is  overcome  by  averaging  the  equations  in  the  limit  T///Ty  — »  0 
wlu're  Tis  is  a  charachaistic  time  describing  eiuTgy  loss  by  convection.  This 
av<Taging  yi<'lds  a  lu'w  theory  from  whi<‘h  a  considerable'  amount  of  infor¬ 
mation  can  b('  de<luc<'d.  Si)ecifically,  the  tenii)erature  is  .si)atially  uniform  in 
the  targi’t  and  ('volv<'s  in  time  according  to  a  first  ord('r,  ordinary  diflerential 
etpiatioii.  The  nonlineanty  in  this  ecpiation  is  a  functional  of  tlu'  ('lei'tric  fic'ld 
within  the  target.  This  equation  is  solved  ftn  a  numlu'r  of  specific  ('xam])les 
and  phy.sical  conclusions  are  drawn  about  certain  heating  processes.  Finally, 
tlu’  problem  of  controlh'd  lu'atiug  is  addre.s.sed  wlu’ie  liiu'ar  feedback  is  shown 
to  be  adequatt'  in  achit'viug  a  jnedeterniiiu'd  tinal  t<'mperature 

4.  Microwa/e  Processing  of  Thin  Ceramic  Rods. 

When  a  thin  cylindrieal  ceramic  sample  is  heated  in  a  single  mode  waveg¬ 
uide  applicator  an  interesting  thermal  phenomenon  ofteji  occurs.  A  locah/;e(i 
hot  spot  forms  near  the  cent<'r  of  the  .sample  and  spreads  outward  from  this 
])oint  raising  the  entire  ceramic  to  an  elevated  tem]>erature  (v-vi).  This  oc¬ 
curs  ('V('!i  though  the  electric  field  intensity  is  es.'^eiil ially  <'onstant  along  the 
axis  of  the  cylinder. 
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We  have  constructed  a  simple  mathematical  [9]  model  based  upon  the 
small  Biot  number  limit  which  gives  a  plausible  explanation  of  the  forma¬ 
tion,  propagation,  and  growth  of  these  hot-spots.  Our  model  is  simpler  and 
easier  to  analyze  than  the  one  used  in  (v)  because  the  variations  of  the 
temj^erature  in  the  cross-section  have  been  averaged  out  and  this  yields  a 
reaction-diffusion  equation  in  only  one  spatial  variable.  The  dimensionless 
diffusion  coefficient  is  given  by 

where  a  is  the  radius  of  the  rod,  L  is  its  length,  and  2?,  is  the  Biot  inunber. 
This  constant  can  be  quite  small  csi)ccially  when  one  is  interested  in  sintering 
ceramic  fibers.  Thus,  the  ju-oldcm  reduces  to  a  singulai’ly  pc'rturbcxl  reaction- 
diffiision  equation  of  which  much  is  known. 

We  have  recently  extended  this  work  along  two  lines.  In  the  first  [14]  we 
have  allowed  the  thermal  conductivity,  the  heiit  capacity,  and  tlu;  density 
of  tlu'  ceramic  to  depend  upon  temperature.  Although  these  dcqumdencies 
are  usually  lu'glected  wlum  compared  to  the  orders  of  magnitude  change' 
that  occurs  in  the  electrical  conductivity,  they  are  .significant  in  affecting 
th('  velocity  and  ultimate  character  of  the  hot-spot  dyni.mics.  A  comi)h't<' 
(h'seription  of  thi.s  mechanism  is  contaiiw'd  in  Reference'  14. 

W('  have  also  analyze'd  this  lu'ating  probh'in  wln'ii  the  sample  is  rotated  90" 
in  tlu'  cavity  [1G|.  The'  mode'l  now  takes  the'  form  of  a  noiiliiu'ar  parabolic 
('(piation  of  reaction-eliffusion  type',  with  a  spatially  varying  re'actiem  te'rm 
that  e'ejrre’spemels  te)  the'  spatial  variatie)n  e)f  the  e'le'e'tre)magne'tic  fie'lel  stre'iigth 
in  the-  wave-guide-.  The-  e-ejuatie)!!  is  analyze-el  aiiel  a  se)luti()n  is  founel  whie-h 
ele'V(-le)i)s  it  he)t  spot  ne-ar  the-  <-e-nte-r  e)f  the'  cylinelrie-iil  s.nnple-  anel  whie-h  the-n 
pre)])agate-s  eiutwarels  until  it  stal)ilize's.  The-  i>ropagatie).->  and  stabilizatie)n 
l)he'nome'ne)n  e-euice-ntrate-s  the-  mie-reaviive- ene-rgy  in  a  le)e-alize'el  re'gie)n  abe)ut 
the'  e-e-nte'r  whe-re-  e-le'vate-el  -mpe-rat ure-s  may  be-  eie-sirable-.  Thus,  we-  eilfe-r  a 
plausible-  e  xi)lauatie)n  fe)r  the-  e-xpe-rime-nts  ele'scribe-el  in  (vii)  whe-re-  lea-alize-el 
fie-lels  iire-  use-el  tee  assist  the-  jeaning  e)f  twe)  e-e-ramie-  re)els. 

Finally,  we-  have-  meiele-le-el  anel  analyze-el  the-  mie  re)wave-  he-ating  e)f  a  e-arbem 
e-e)ate-el  e-e’remiie-  fil)e-r  [15,17].  Sine-e-  the-  e'le-e-trie-al  e-e)ndue  tivity  e)f  the-  e-arbem 
is  nme-h  gre-ale-i  than  that  e)f  the-  ce'ramie-  at  re)e)m  -ni]>e-r!iture',  the-  e-e)iiting 
rapidly  he-ats  einel  transfe-rs  the-ruuii  e-ne-rgy  to  tl  e-e-ramie-  ley  e-eenelue-tieen. 
This  raise-s  tlie-  ave-rage-  te-mpe-rat are'  of  the-  mate-rial,  een  the-  e-eenve-e-tive-  time- 
se-ale-,  anel  affe-e-ts  the-  initial  e-eenehtieen  e)f  eatr  small  Bieet  the-eery.  The-  re-sults 
pie-elie-t e-el  by  the'  the-eery  ejualitative-ly  agre-«-  with  e-xpe-rime-nl s  (viii)  anel  give- 
a  ue-w  ine-e  hiuii.-an  fe  r  re-aehing  hithe-rle)  inae-e-e-ssilele-  leeertieens  eif  the-  s-shaiee-el 
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response  curve.  A  complete  discriptioii  of  this  mechanism  is  contained  in 
Reference  17. 

5.  Asymptotic/Nuinerical  Hybrid  Methods:  An  Example  from 
Acoustic  Scattering  Theory. 

We  have  developed  a  hybrid  method  [3],  based  upon  the  method  of  matched 
asymptotic,  expansions,  to  study  the  scattering  of  a  plane  acoustic  wave  from 
a  large  finite  diffraction  grating.  The  inner  problem  of  an  infinite  grating 
is  solved  by  a  well-posed  numerical  method  from  which  the  modal  ampli¬ 
tudes  are  deternrincd.  The  outgoing  scattered  plane  waves  do  not  satisfy  the 
Somerfeld  condition  so  that  this  approximation  becomes  nonuniform  in  the 
far  field.  The  outer  problem  of  a  finite  strij)  of  unknown  sources  is  approx¬ 
imated  by  the  method  of  stationary  ph:isc  and  the  results  are  matched  to 
those  of  th(?  inner  problem.  The  matching  determines  the  amplitudes  of  the 
unknown  sources.  We  have  i)rcsented  an  explicit  formula  for  the  scattering 
cross  section  which  blends  Ixith  numerical  and  analytical  results.  It  shows 
a  highly  localised  scattering  pattern  whose  maxima  coincide  with  the  Bragg 
angle's  of  the  infinite'  perie)dic  structure.  This  localization  as  well  as  the  niag- 
nituele's  of  the'  maxima  e'le'arly  eleinoustrate  the'  effects  e)f  the  grating’s  finite 
size'. 

6.  Large  Finite  Dielectric  Gratings. 

The'  i)ioble'm  e)f  plane'  wave  scattering  by  a  finite',  but  ('lectrically  large, 
l)erie)die'  elie'h'ctrie'  grating  on  the-  inte'rfae'e  be'twe'cn  se'mi-infinite  dielectric 
me'elia  is  an  imi)e)rtant  nie)ele'l  of  a  microwave' e'lU'rgy  e'o)ipl('r.  When  the'  grat¬ 
ing  parame'te'is  fall  in  the'  re'.seaimie'e'  re'gime',  i.e'.,  the'  he'ight  and  pe'riexl  of 
the-  grating  are'  comparable'  te)  the'  wave-length  eef  the-  ine'ieie'ut  wave,  them  nu- 
me’rie'al  me'theeels  l»e'ce)nu'  ina])pre)priate'  be'e'au.se'  eef  the'  eli.si)arate'  h’Ugth  scale's 
l)re'.s<'nt.  We  have  ajjplie'el  die'  Me-theiel  e)f  .Matche-el  Asymi)te)t ie'  Expansiems 
tei  e-emstrue-t  a  hyliriel  nunu'rie'al/a.syinptotir  me'theiel  feir  e-ificie-utly  analyzing 
this  preible'Ul  [IS]. 

.AiJiiroiu'iate'  scalings  eif  the-  inele-pe-nele-nt  varialile's  liavi-  re-sulte-el  in  iuiie-r 
anel  eiute-r  preible-ms.  wlieise-  re-spe-e-tive-  .solufiems  are'  eibtailie'ei  liy  a  luune'ri- 
e-al  and  an  analytiesil  me-thoel,  re-spe-e-Mve-ly.  The-  appreiaeh  has  proviele'd  an 
alte'inative  mi'theiel  to  le-ngthy  numerie  al  e-ale'ulatieais  anel  has  maele-  the-  e’har- 
ae'te'i'izat iem  eif  se-atte  riug  liy  finite'  grating.-  whei.se-  li-ngth  is  of  the-  eireler  eif 
te-ns  of  wavi'le'iigt hs.  A  eleiseil  form  e-xpressieiji  for  the-  scatte'iing  patte-rn 
has  be-e-n  lie-lived  in  ti-rnis  eif  the-  Hayle-ig'.i  re-flection  anil  transmission  moehil 
ainplit  lull's  whie'h  have-  bi-i'ii  olitaine-el  from  the-  nuiMi'rii-al  approai'h  u.se-el  to 
solve  the'  inner  jirolili'm.  Tin-  ele-tails  of  the-.se'  I'alculations  aiul  ri'sults  are' 
containe'il  in  He'fe'ie'ni’e'  18. 


7.  Large  Membrane  Array  Scattering. 

Wc  have  extended  the  above  techniques  to  study  the  scattering  of  a  plane 
acoustic  wave  from  a  large  array  of  baffled  membranes  [10].  The  structure 
is  a  simple  model  of  an  acoustic  antenna.  The  far  field  pattern  again  is 
highly  localized  aljout  the  Bragg  or  mode  angles  of  the  inner  representation. 
The  maximum  value  of  this  function  is  directly  proportional  to  the  reflection 
coefficient  determined  numerically  from  the  inner  problem.  Several  numerical 
examples  are  presented  illustrating  these  features.  Moreover,  it  is  observed 
that  the  asymptotic  theory  agrees  qxiite  well  with  the  exact  answer  (obtained 
by  using  a  finite  difference  scheme)  when  there  are  only  three  membranes 
present . 

Wc  have  also  developed  a  new  numerical  method  for  solving  the  inner  prob¬ 
lem  (infinite  grating)  when  many  Bragg  modes  are  present.  The  method  is 
based  upon  a  finite  difference  scheme  used  in  conjunction  with  a  new  radi¬ 
ation  boundary  operator.  This  operator  is  local  in  character  and  effectively 
allows  the  Bragg  modes  to  radiate  out  of  the  computational  domain  (funda¬ 
mental  cell).  A  full  description  of  the  method  is  given  in  Reference  13. 

8.  Acoustic  Target  Reconstruction  Using  Geometrical  Optics  Phase 
Information. 

A  simple  algorithm,  which  uses  the  phase;  information  from  the  geomet¬ 
rical  optics  limit,  to  construct  tlie  shape  of  an  object  has  been  presented 
in  R(;forcncc  5.  Essentially,  the  phase  in  the  back  scattered  direction  was 
shown  to  determine  the  equation  of  the  tangent  plane  at  the  unique,  but 
unknown,  specular  point.  This  plane  depends  upon  the  two  spherical  angles 
9  and  <jf>,  which  describe  the  incident  wave  direction.  The  observation  that 
tin*  tang<;nt  plaiu;  envelopes  the  obstacle  as  9  ami  <t)  are  varied,  allowed  the 
derivation  of  an  explicit  formula  for  the  equation  of  the  surface  in  terms  of 
the  measured  .sca.ttor<;d  phas('.  An  analogous  two-dimensional  formula  has 
also  been  prc'sented.  A  full  description  of  the  nu'thod  is  given  in  Rehn'cnce 
5. 

9.  Rapid  Pulse  Responses  for  Scattering  Problems. 

We  have  recently  considered  the  numerical  computation  of  the  fi{;ld  scat¬ 
tered  from  a  body  in  two  dimensions  due  to  an  incident  plane  pressure  pulse. 
In  particular,  W(;  have;  examined  the  process  of  inferring  the  scattered  field 
du<'  to  ou('  incident  pulse  given  tin'  scattered  field  due  to  another  incident 
pulse'.  Tin'  oljjective  is  to  devt'lop  an  indirect  method  that  avoids  the  poten¬ 
tially  ('xp('nsiv('  direct  solution  of  the  problem.  Our  apinoach  is  based  on  a 
fornuda  c'xpn'ssing  tin'  scattc'ii'd  field  as  a  <-onvolution  of  a  kc'rind  with  the  in¬ 
cident  puls('  profile'.  The  most  straight  forward  gein'ralization  of  this  formula 


S 


to  the  discrete  version  of  tlie  scatter  problem  used  iii  numerical  computations 
does  not  allow  the  k(;rnel  to  be  inferred  from  a  single  numerical  experiment — 
a  difficulty  we  have  called  the  multi-source  problem.  Preprocessing  the  inci¬ 
dent  pidses  using  simple  interpolation  formuhis  overcomes  the  inuUi-soui'ce 
problem  giving  an  exact  algorithm  for  computing  the  kernel.  Selection  of 
a  sharp  incident  pulse  (the  Kronecker  pulse)  for  the  primary  numerical  ex¬ 
periment  assures  stability  of  this  algmithm  and  permits  extremely  accurate 
prediction  of  the  scattered  fields  for  secondary  incident  pulses.  A  complete 
description  of  our  results  can  be  found  in  Reference  11. 
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